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ABSTRACT 

We calculate a new equation of state for baryons at sub-nuclear densities 
meant for the use in core-collapse simulations of massive stars. The abundance 
of various nuclei is obtained together with the thermodynamic quantities. The 
o ,| formulation is the NSE description and the liquid drop approximation of nuclei. 

q The model free energy to minimize is calculated by relativistic mean field the- 

ory for nucleons and the mass formula for nuclei with the atomic number up to 
~ 1000. We have also taken into account the pasta phase, thanks to which the 
transition to uniform nuclear matter in our EOS occurs in the conventional man- 
ner: nuclei are not dissociated to nucleons but survive right up to the transition 
O's ■ to uniform nuclear matter. We find that the free energy and other thermodynam- 

i—{ ; ical quantities are not very different from those given in the H. Shen's EOS, one 

of the standard EOS's that adopt the single nucleus approximation. The average 
mass is systematically different, on the other hand, which may have an important 
t— i ■ ramification to the rates of electron captures and coherent neutrino scatterings 

on nuclei in supernova cores. It is also interesting that the root mean square 
of the mass number is not very different from the average mass number, since 
the former is important for the evaluation of coherent scattering rates on nuclei 
but has been unavailable so far. The EOS table is currently under construction, 
which will include the weak interaction rates. 
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Introduction 



Equations of state (EOS's) of hot and dense matter play an important role in the 
dynamics of core collapse supernovae both at sub- and supra-nuclear densities (see e.g. 
Janka et al.l (120071 ) ) . Not only thermodynamical quantities such as pressure, internal energy, 
entropy, sound velocity but also information on matter composition is provided by EOS. At 
sub-nuclear densities, the latter affects the rates of electron captures and coherent neutrino 
scatterings on nuclei, both of which in turn determine the electron fraction of collapsing 
cores, one of the most critical ingredients for the core dynamics. This is exactly the issue we 
would like to address in this paper. 

The EOS employed for the core collapse simulation must cover a wide range of density 
(10 5 < pb < 10 15 g/cm 3 ) and temperature (10 9 < T < 10 12 K), including both neutron- 
rich and proton-rich matter. One of the difficulties in constructing the EOS originates 
from the fact that depending on the density, temperature and proton fraction, the matter 
consists of either dilute free nucleons or a mixture of an ensemble of nuclei and free nucleons 
or strongly interacting dense nucleons. Another complication is the existence of the so- 
called nuclear pasta phase, in which nuclear shapes change from droplet to rod to slab to 
anti-rod and bubble(anti-droplet) as the density i ncreases toward the nuclear saturation 
density, at which uniform n uclear matter is realized ( iRavenhall et al.lll983t lHashimoto et al. 
1984 ; iNakazato et al.l 120091 ) . At high temperatures (T>5x 10 9 K), chemical equilibrium is 
achieved for all strong and electromagnetic reactions, which is referred to as nuclear statistical 
equilibrium, or NSE, and the nuclear compos ition in the matter is determined as a function 
of density, temperature and proton fraction. (IBlinnikov et al.ll2009l ) At lower temperatures, 
the matter composition is an outcome of preceding nuclear burnings and cannot be obtained 
by statistical mechanics. In this paper we are concerned with the high temperature regime, 
in which the nuclear composition is a part of EOS. 

At present, there ar e only two EOS's in wide use for the core-collapse simulation. The 
Lattimer-Swesty's EOS (ILattimer et al.lll99ll ) is based on Skyrme-type nuclear interactions 
and the so-called compressible liquid drop model for nuclei surrounded by dripp ed nucle- 
ons. The EOS by H. Shen, Toki, Oyamatsu and Sumiyoshi (IH. Shen et al.lll998l ) employs 
a relativistic mean field theory (RMF) to describe nuclear matter and the Thomas-Fermi 
approximation for finite nuclei with dripped nucleons. It should be emphasized here that 
both EOS's adopt the so-called single nucleus approximation (SNA), in which only a sin- 
gle representat i ve nu cleus is included. In other words, the distribution of nuclei is ignored. 
Burrows et al.l (119841 ) demonstrated that SNA is not a bad approximation for thermodynam- 
ical quantities such as pressure. It will not be the case, however, for the weak interaction 
rates, since the electron capture rates are sensitive to nuclear shell structures and the great- 
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est contributor is not the mos t abundant nuclei tha t the single representative nuclei in SNA 
are supposed to approximate (ILanganke et al.ll2003l ). It is also noted that the coherent scat- 
tering of neutrinos on nuclei has a cross section proportional to A 2 with A being the nuclear 
mass number. In SNA we have no information on the average of mass number squared A 2 
and replace it by A 2 , with A rep being the mass number of the representative nucleus. This 
is certainly an approximation that needs justification. 



Recently iHempel et al.l (120101 ) published a new EOS at sub- nuclear densities meant 
for the use in supernova simulations. Assuming NSE, they calculated the abundances of 
nuclei as a function of density, temperature and proton fraction, minimizing the model free 
energy. They employed RMF for the nucleon contribution and adopted masses of nuclei 
from two mass tables in pu blic domain. One of them is the experimental atomic mass 
table 2003 flAudi et al.l 120031 ) and the other is a tabulated theo retical prediction for nuclei 
with no available exp erimental data, which is based on RMF (IGeng et al.l 120051 ) with the 
TMA parameter set (IToki et al.lll995l ). The translational and Coulomb energies of nuclei 
are designed so that the free energy could match the one given by RMF for uniform nuclear 
matter at the saturation density. It should be not ed that the solutio n of the Saha equation for 
the ma ss data, which are obtained in laboratory (lAudi et al.ll2003l ) or calculated for isolated 
nuclei (IGeng et al.ll2005l ). is not appropriate at high densities (ps > 10 12 g/cm 3 ), since nuclei 
are not isolated and affected by surrounding other nuclei, nucleons and electrons, resulting 
in the modification of masses of nuclei. 



Our formulation is quite similar to the one employed by lHempel et al.l (|2010| ). The major 
difference, however, is that we calcula te very heavy nuclei (Z > 100 with Z being the atomic 
number) with a mass formula whereas IHempel et al.l (120101 ) ignored these nuclei owing to the 
limitation of the ma ss tables. It is stress ed, however, that these massive nuclei do appear in 
the H. Shen's EOS (|H. Shen et al.lll998l ) and, as shown later, they are indeed abundant at 
high densities (ps > 10 13 g/cm 3 ) in our EOS also. Another important improvement in this 
paper i s the treatment of t he nuclear pasta phase, which was entirely ignored in the previous 
paper (IHempel et al.ll2010l ). It turns out that this is important to reproduce the conventional 
manner of transition to uniform nuclear matter, that is, nuclei do not melt down to nucleons 
but survive up to the nuclear saturation density. 

The aim of our project is in a sense to merge all the previous EOS's to provide all the 
information needed for core collapse simulations: thermodynamical quantities such as pres- 
sure, internal energy, entropy, sound velocity and some other important derivatives as well as 
information on the abundance of nuclei required for the calculation of the weak interaction 
rates; the description of uniform nuclear matter is based on RMF in H. Shen's EOS; the ap- 
proximate treatment of the pasta phase is similar to the one employed in Lattimer-Swesty's 
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EOS; we solve a Saha-like equation to obtain the abundance of nuclei. In this paper, we 
describe the formulation in detail and show some features of our EOS, which are particularly 
different from those of the previous EOS's. More systematic and thorough presentation of 
the EOS will be given in the sequel after we have constructed the EOS table that includes 
somehow the electron capture rates. 

This article is organized as follows. In section 2 we explain in detail the model free 
energy to be minimized. Then the results of minimization are shown in section 3, with an 
emphasis being put on the comparison with the previous EOS's. The paper is wrapped up 
with a summary and some discussions in section 4. 

2. Model Free Energy 

The strategy to obtain the multi-component EOS is to construct a model free energy and 
minimize it with respect to the parameters involved in the expression of the free energy. The 
matter in the supernova core at sub-nuclear densities consists of nucleons, nuclei together 
with electrons and photons. The latter two are not treated in this paper although the 
inclusion of them as ideal Fermi and Bose gases respectively are quite simple and now a 
routine. Neutrinos are not always in thermal or chemical equilibrium with the matter and 
cannot be included in the EOS. Their non-equilibrium distributions should be computed by 
the transport equations. 

The free energy is thus a sum of the contributions from the individual nucleus and those 
from the nucleons that are not contained in the nuclei. We employ the RMF free energy 
for the nucleons outside the nuclei. In constructing the model free energy for the nuclei, the 
following points are appropriately taken into account: 

1. Not only the nuclei that are known to exist in laboratory but exotic ones such as 
very neutron-rich and/or very heavy nuclei should be also included somehow even if 
experimental data are unavailable. They are indeed obtained in supernova simulations 
with EOS's that employ SNA. 

2. At low densities, experimental mass data should be used whenever available. This 
is important to guarantee that the most stable iron group nuclei are populated with 
correct abundances at low temperatures. 

3. At higher densities, strong and electromagnetic interactions among free nucleons, nuclei 
and electrons cannot be ignored and should be properly reflected in the surface and 
Coulomb energies of nuclei. 
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4. The bulk energy of nuclei should approach the one for uniform nuclear matter as 
the density increases to the nuclear saturation density. This should be taken into 
consideration to obtain a continuous transition to uniform nuclear matter. 

5. Near the nuclear saturation density, the nuclear pasta phase should be treated some- 
how. This is also a necessary ingredient to realize a seamless transition to uniform 
nuclear matter, in which nuclei are not dissociated to nucleons but continue to exist 
right up to the transition. 

The model free energy for nuclei described in detail in the following subsections is 
based on the mass formula. As mentioned above in point [2J at low densities, experimental 
mass data are utilized, whenever available, to obtain the bulk energy of nuclei and RMF is 
employed otherwise. At high densities, however, even for the nuclei, for which experimental 
data exist, the bulk energy is modified so that points [3] and @] are satisfied. As for point El 
only the bu bble phase is explicitly considered, following the phenomenological prescription 



adopted in (lLattimer et al.l Il99fl ). Although it is easy to change it, we suppose in this 
paper that H. Shen's EOS is used for uniform nuclear matter at supra-nuclear densities. 
The continuous transition from the EOS at sub-nuclear densities to the one at supra-nuclear 
densities is guaranteed by the following prescriptions: the free energy of the nucleons outside 
nuclei is based on the same RMF model as for uniform nuclear matter; the bulk energy of 
nuclei is modified to approach that of uniform nuclear matter as the density increases; the 
surface, Coulomb and translational energies of nuclei go to zero as the transition density is 
approached. 

The details of these prescriptions will be given in the following subsections. We briefly 
explain the RMF model we employ for the nucleons in section 12.11 Then we proceed to the 
description of the model free energy for nuclei in section I2.2[ in which the Coulomb and 
surface energies are treated in subsections 12.2.11 and 12. 2. 2} respectively, the bulk energy is 
presented in subsection 12.2.31 and the translational energy is given in subsection 12.2.41 The 
method to minimize the model free energy is explained in section 12.31 and the evaluation of 
thermodynamical quantities from the free energy is mentioned in section 12.41 



2.1. Free Energy for Nucleons outside Nuclei 



Th e free energy of the nucleons that are not contained in nuclei is base on RMF in this 
paper (|H. Shen et al.lll998l ). The choice is not essential, however, and one can replace it by 
any other theory or model indeed. The important thing to ensure the continuous transition 
to uniform nuclear matter is to use the same theory or model for the nucleons outside nuclei 
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at sub-nuclear densities and uniform nuclear matter at supra-nuclear densiti es. The details 



of the RMF model employed in this paper can be found in the original paper (lH. Shen et al. 
19981 ). In the following, we summarize only the essential ingredients. 



In RMF nuclear interactions are described by the exchange of mesons. The employed 
Lagrangian is the following: 



£rmf = - M - g a a 
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In the above expression, the notation is the same as in |H. Shen et al.l (119981 ): ip, a, u and 
p denote nucleons (proton and neutron), scalar-isoscalar meson, vector-isoscalar meson and 
vector-isovector meson, respectively, and W pv = d^uj u — d u u tl and K a = d fl p ai ' — d u p atl + 
The nucleon-meson interactions are given by the Yukawa couplings and the 
isoscalar mesons (a and u) interact with themselves, which are expressed as the cubic and 
quartic terms. In the mean field theory, the mesons are assumed to be classical and re- 
placed by their ensemble averages whereas the Dirac equation for nucleons is quantized and 
the free energy is evaluated based on the energy spectrum of nucleons obtained this way. 
M is the mass of nucleon s and assumed to be 938MeV. We use the TM1 parameter set 



as m 



H. Shen et al.l (119981 ). in which the masses of mesons, 



ni 



and the coupling 

constants, g a , g W) g p , g 2 , #3, c 3 , are determined so that not only the saturatio n of uniform nu 



clear matter but also the properties of finite nuclei could be best reproduced (ISugahara et al. 
I99J). 



It is interesting to study how the resulting EOS is affected by the choice of RMF param- 
eter or , more generally, by the choice of EOS for uniform nuclear matter. If the preceding 



paper (IHempel et al.ll2010f ). in which several RMF parametrizations were compared, is any 
guide, differences in the RMF parametrization will be reflected in thermodynamic quantities 
near the saturation density. Since RMF is employed also for the bulk energy of nuclei in our 
EOS as explained later in subsection 12.2.31 differences in symmetry energy originated from 
different parametrizations may be reflected in the neutron-richness of nuclei, for which we 
have no available experimental mass data and the empirical mass formula is employed. The 
abundance of nuclei with experimental mass data, on the other hand, will not be changed very 



much. It should be also noted that ! 
ter energy at sub-nuclear densities ( 


-IMF in general 


;ends to overestimate the neutron mat- 


Sumivoshi et al. 


1995; 


Akmal et al. 


1998; 


Steiner et al. 


2005; 


Gandolfi et al. 


2010; 


Hebeler et al. 


2010 


). It is stressed that one of the features of our 
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EOS is that it is easy to adopt different EOS's at supra-nuclear densities as long as they are 
available. The above issues will be discussed more in detail in the forthcoming paper. 

In this paper, the free energy density of uniform nuclear matter, f RMF (riB,T,Y p ) with 
ns,T and Y p being the baryon number density, temperature and proton fraction, respectively, 
is used mainly at sub-nuclear densities. In so doing, one important thing to be taken into 
account is the fact that the free nucleons cannot exist in the volume occupied by nuclei. The 
volume, V, that can accommodate the nucleons is hence given by 

V^V-Y^NV", (2) 

i 

where index i specifies a nucleus, V is the total volume, iVj is the particle number of nucleus 
i and is the volume of nucleus i, which is expressed by the mass number, Ai, and number 
density, rii, of the same nucleus as well as the nuclear saturation density, n S i, as follows: 

Vf = (3) 

where we assume that the nuclei are homogeneous, having the nuclear saturation density, 
n si . The reason why n si has the subscript i will be made clear later (see section |2~2|) . Then 
the densities of free protons and neutrons in the volume V, which are denoted by n p and 
n' n , respectively, are given by 

where N p / n and n p / n are the particle numbers and the number densities in the total volume 
V of nucleons and rj is expressed as rj = V /V = 1 — Y2i^-i n i/ n si- The free energy density 
for the nucleons in V should be evaluated at Y p = n' p /{n' p + n' n ) and n B = {n' p + n' n ) and is 
given for the entire volume as 

fp,n = V fp,n (5) 

withf; :n = f RMF (n' p ,n' nJ T). 



n P /n — ^rr — W 



2.2. Free Energy for Nuclei 

This is the most crucial part of our model free energy and is described in detail in 
the following. The free energy of nuclei consists of two parts, one of which concerns a 
translational energy whereas the other part is roughly a nuclear mass. The modeling of the 
latter is based on a semi-empirical mass formula and takes into account the Coulomb, surface 
and bulk energies. We explain each contribution to the mass formula one by one and then 
move on to the treatment of the translational energy. 
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Before proceeding, however, we state here the definition of the nuclear saturation density, 
n S i(Zi/Ai, T), more precisely, since it will appear frequently in the following. It is in principle 
the baryon number density at which the free energy per baryon, F RMF (T,nB,Y p ) with 
Y p = Zi/Ai, by RMF takes its minimum value or, put another way, the uniform nuclear 
matter with the same proton fraction as the nucleus has the lowest energy in RMF. Note 
that the proton fraction of each nuclei is different from the one for the whole system and 
that our nuclear saturation density is actually different from nucleus to nucleus and hence 
has the subscript that indicates which nucleus is referred to. There are two exceptions to 
this definition, though. The first one is the case, in which the saturation density given above 
is smaller than the number density of the whole system, n#. This happens if ns reaches 
the saturation density of some nucleus, n S i(Zi/Ai,T), before the whole system becomes 
uniform. This is indeed possible because the saturation density tends to be smaller for more 
neutron-rich nuclei. In such a case the nucleus will be overcompressed and we replace the 
original saturation density with the number density of the whole system, ns- The second 
case occurs at high temperatures, for which the free energy, F RMF (T,nB, Zi/Ai), has no 
minimum because the entropy contribution, the TS term with S being entropy, overwhelms 
the internal energy. In that case the nuclear saturation density, n si , is set to be the density 
where the matter becomes uniform in the H. Shen's EOS. 



2.2.1. Coulomb energy 

To calculate the Coulomb energy density in the free energy for nuclei we introduce a 
Wigner-Seitz cell (W-S cell) for each species of nuclei and impose charge neutrality in it. It 
is assumed that a nucleus is centered in the W-S cell, the volume of which is denoted as Vi 
and that the cell contains also nucleons as a vapor outside the nucleus as well as electrons, 
which are assumed to be uniform in the entire cell. The charge neutrality in the cell can be 
expressed as 

V<n e = Z t + (Vi - V?)n' p , (6) 
which can be solved for Vi and gives 

n e - n' p 

where Vf is the nuclear volume in the cell and n e denotes the number density of electrons. 
The vapor volume and nuclear volume fraction for each nucleus are given by V R = Vi — V^ 
and u,i = V { N /Vi, respectively. 

As the density approaches the nuclear saturation density, the nuclear pasta phase is 
encountered, in which various shapes of nuclei exist. As mentioned already, this phase has 
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to be somehow treated in order to obtain the smooth transition to uniform nuclear matter. 
It is much beyond the scope of our project to treat all known pastas accurately, however. 
We hence resort to a rather crude approximation, in which the bubble phase is explicitly 
considered. It is known that bubble nuclei, in which the vapor nucleons are surrounded by 
the nucleus, appear at the end of the pasta phase. We assume there that a spherical vapor 
is located at the center of W-S cell surrounded by a shell-shaped nucleus. In this paper 
we assume further that each nucleus enters the nuclear pasta phase individually when the 
yolume fraction, iij, re aches 0.3 and that the bubble shape is realized when it exceeds 0.7 
( IWatanabe et al.l 120051 ). The intermediate states (0.3 < U{ < 0.7) are simply interpolated 
from the normal and bubble states. This criterion is admittedly rather arbitrary. We have 
tried another choice, 0.4 < Ui < 0.6, and confirmed that the results are hardly affected, 
though. We hence believe that our EOS is not very sensitive to the criterion. 

The evaluation of the Coulomb energy in the W-S cell is straightforward and given by 
the following expression: 

E° = [ g(r)dg(r) , (8) 

J cell r 

where the cell is assumed to be spherically symmetric and q(r) is the sum of electron and 
proton charges up to the radius r. Assuming further that the protons are uniformly dis- 
tributed in the nucleus and the vapor separately, we obtain the Coulomb energies for the 
normal and bubble phases from the above integration as follows: 
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{ui > 0.7), 



(9) 



with D(ui) = 1 — |m^ 3 + \ui, where e is the elementary charge. Note that the two energies 
coincide with each other at Ui = 0.5, since the nuclear volume, V/^ , is equal to the vapor 
volume, Vj 8 , there although the value is not used and the interpolation is done from the 
values at Uj = 0.3 and Ui = 0.7 to take crudely into account the existence of other pastas. 
We employ cubic polynomials of Ui for interpolation. The four coefficients of the polynomial 
are determined by the condition that the Coulomb energy is continuous and smooth as a 
function of at m = 0.3 and u,- = 0.7. 
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2.2.2. surface energy 



In the semi-empirical mass formula we employ in this paper, the surface energy of 
nuclei is given by the product of the nuclear surface area and the surface tension, which is 
approximated by the following analytic expression: 



.4 



2/3 



Pi 



O"0 



Airr 2 N . 



[S s (l-2{Z i /A i ] 



(10) 



where r Ni = (3/47rV^ Ar ) 1 / 3 is the radius of nucleus i and a denotes the surface tension for 
symmetric nuclei. The second term on the right hand side is added to accommodate highly 
neutron-rich nuclei, for which no experimental data are available but a proper consideration 
of the surface symmetry energy is almost mandatory. The values of the constant s , an = 
1.15MeV/fm 3 and S s = 45.8MeV, are adopted from the paper by lLattimer et al.l ( Il99ll ). 



Note that this expression gives negative values if nuclei become extremely neutron-rich. We 
hence assume in this paper that only those neutron-rich nuclei with a surface tension higher 
than a prescribed lower limit, o = 0.272MeV/fm 3 , can exist. This yalue is obtained from 
the mass data of the most neutron-rich nuclei with {Z,A) = (1,7) ( lAudi et al.l 120031 ) and 
the subtraction of the bulk and Coulomb energies given by our mass formula. 

Noting that the nuclear surface should be replaced by the bubble surface in the bubble 
phase, we obtain the following expression for the surface energy: 



E 



4irr 2 m Oi I 1 v - 



1-/7,-, n, ( 1 



n'+n 



I \ 2 



2/3 



(Till- 



2/3 

Itt ( ^-V? ) a, f I 

4:71 



n' + n 



I \ 2 



p 1 n 



n 'p + n n 



I \ 2 



fa < 0.3), 
(ut > 0.7), 



11) 



where r Bi = (S/AttV^) 1 ^ 3 is the radius of bubble i. The last factor, (l — (n' p + n' n )/n si ) , in 
these expressions are introduced to take into account the effect that the surface energy should 
be reduced as the density contrast decreases between the nucleus and the nucleon vapor. The 
specific functional form is inspired by the surface energy empl oyed in the T h omas -Fermi 
approximation in the H. Shen's EOS, |V(n n + n p )\ 2 (Eq. (35) of |H. Shen et al.l ( 19981 ) ). As 
shown later, this prescription is necessary to ensure the continuous transition to uniform 
nuclear matter, since the transition occurs not only by coalescence of nuclei but also by the 
increase in the density of nucleon vapor. For the other pastas (0.3 < Ui < 0.7), we employ 
the same cubic interpolation as for the Coulomb energy (section I2.2.ip . 
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2.2.3. bulk energy 

The bulk energy of nuclei should be different at high densities, ps > 10 12 g/cm 3 , from 
the one at lower densities, since nuclei cannot be regarded as an isolated ob ject and effects 



of other nuclei, nucleons and electrons cannot be ignored (IVautherin Ill994l ). In this paper, 
we will take this into account phenomenologically as follows. 

As the density increases, the nuclear bulk energy, Ef, should approach the one for 
uniform nuclear matter with the same proton fraction. We obtain the latter from RMF as 

Ef{T) = A t F RMF (n sl , T, Zi/Ai), (12) 

where F RMF \ub, T, Y p ) is the free energy per baryon by RMF. Note that the proton fraction 
in this expression is not the one for the whole system but the one for each nucleus; the 
density is taken to be the saturation density defined earlier (section 12. 2p for each nucleus. 
This implies that nuclei have their own densities that are different from each other. Since the 
saturation density tends to be lower for neutron-rich nuclei, the density of the whole system 
becomes equal to the one for some neutron-rich nucleus one after another as it increases. In 
such a case the density of the nucleus that becomes uniform in its W-S cell is set to be the 
system density as mentioned earlier. These prescriptions ensure the continuous transition to 
uniform nuclear matter. 

At low densities, o n the other hand, experimental mass data are used whenever avail- 



able (I Audi et al.l 120031 ). This is indeed important to reproduce the well-known results of 
the ordinary NSE. As will be shown later, the nuclear shell and even-odd effects manifest 
themselves in the abundance of nuclei and should be taken into account whenever possible. 
These experimental mass data are accommodated in our mass formula by adjusting the bulk 
energy, which is expressed as 

E^ 1 vacuum ] vacuum i (l*^) 

where subtracted from the experimental mass, M^ ata , are the Coulomb and surface energies 

c s 

[E i ' } V acuum of the nucleus isolated in vacuum, that is, the energies calculated from our mass 
formula for n e = n' p = n' n = 0. This ensures that the experimental mass is reproduced 
exactly in the low density limit. At finite densities, on the other hand, the mass is given by 

Mi = Ef + Ef + Ef 

= Mf ata + E?(n! n , n' p , n e ) + E?(n' p , n e ) - [E? } mcuum - [Ef} vacuum . (14) 



Note that the bulk energy in our mass formula actually includes other energies such as the 
bulk symmetry, shell and pair energies. 
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For nuclei with no available experimental data, we employ the bulk energy obtained 
from RMF even at low densities as explained above. Since these nuclei, which are either 
very heavy Zj > 100 or highly neutron-rich Zi/Ai < 0.35, are abundant normally at high 
densities (p# > 10 12 g/cm 3 ) alone, the above prescription is justified. The only exception is 
the combination of very low temperatures (T < lMeV) and low proton fractions (Y p < 0.3). 
Such cases are not encountered in the simulations of core-collapse supernovae, however, and 
will not be a problem practically, either. These low- and high-density values (Eqs. (TLB"]) 
and ( II 2p ) are linearly interpolated between ps = 10 12 g/cm 3 and the nuclear saturation 
density, ps ~ 10 14,2 g/cm 3 . The linear interpolation makes the free energy not smooth 
and the pressure discontinuous at the boundaries of the interpolation region. In practice, 
however, the variation of the bulk energy is quite minor compared with those of Coulomb 
and translational energies and the discontinuities of the pressure are negligible. Moreover, 
nuclei that have experimental mass data and hence need the interpolation are not abundant 
near the high density end, the fact which makes the discontinuities even less important there. 
We hence believe that the above prescription is sufficient. 



2.2.4- translational energy 

The translational energy of nucleus i in our model free energy is based on that for the 
ideal Boltzmann gas and given by 

El=k B TLJ^-)-l)(l-^), (15) 



9i(T)n Qi J J V n 

where fee is the Boltzmann constant and uqi = (MikBT/2Trh 2 ) 3 ^ 2 . Since nuclei have internal 
degrees of freedom, the contribution from excited states of nuclei to the free energy cannot 
be ignored at high tempe ratures. In the a bove formula, the effect is encapsulated in gi(T), 



which is approximated by lFai et al.l (119821 ) as 



9i (T) = 0° + J75 J™ dEe- E / T e W (^2aJX0) , (16) 

where g® denotes the contribution of the ground-state (the spin degrees of freedom and 
normally much smaller than g.i{T)) and is set to be g® = 1 for even nuclei and g® = 3 for 
odd ones in this paper for simplicity. Other coefficients in the contribution from the excited 
states are given as a(A { ) = (A/8)(l - c 2 4~ 1/3 )MeV-\ c x = O^MeV" 1 and c 2 = 0.8. In 
order to estimate the sensitivity of results to gi(T), we multiplied the original g%{T) rather 
arbitrarily by a factor of 0.5 or 2 at Y p = 0.3 and T = 5MeV and compared the results. 
We found small differences at p ~ 10 13 g/cm 3 . For example, the free energy is decreased or 
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increased by 4% at most; the entropy changes -1% and +3% whereas the mass fraction of 
all nuclei varies -10% and +8% from the original values. 

The last factor on the right hand side of Eq. (|T5|) takes account of the excluded-volume 
effect: each nucleus can move in the space that is not occupied by other nuclei and free 
nucleons. The factor reduces the translational energy at high densities and is important to 
ensure the continuous transition to uniform nuclear matter as mentioned earlier. In reality, 
however, the nuclear translational energy is supposed to be suppressed by the formation of 
Coulomb lattice when the so-called plasma factor, which is defined as T = (Ze) 2 /(fk B T) 
with th e average proton nu mber Z and distance between nuclei f, reaches a critical value 



~ 171 (ISlattery et al.l Il998l ) and the introduced factor may be regarded as a very crude 
phenomenological approximation to this situation. The present form of the factor, (1 — 
n B/ n s) — (V — Vbaryon)/V, actually gives a linear suppression in terms of the occupied 
volume Vbaryon- Here we always employ the nuclear saturation density for symmetric nuclei 
n s = [n si (Zi/Ai,T)} Zi / Ai=05 for numerical convenience, which will be justified, provided 
the approximation itself is very crude. Note also that the translational energy is a minor 
contribution except at low densities, where the factor is almost unity. 



2.3. Minimization of Free Energy 

The abundances of nuclei as a function of p B , T and Y p are obtained by minimizing the 
model free energy derived so far with respect to the number densities of nuclei and nucleons 
under the constraints, 



+ n n + ^2Aini = n B =p B /m B , 

% 

+ ^ Z i U i = n e = Y p n B , (17) 



just as in the ordinary NSE calculations. Note that n p / n is the number density of pro- 
ton/neutron in the total volume V of the system as defined in 12. 11 By introducing Lagrange 
multipliers, a and /3, for these constraints, the minimization condition is given by 

d 

{/ - a{n p + n n + y^ j AjUj - n B ) - P{n p + ^ ZjUj - n e )} = 0, (18) 



any- 



where/ is the free energy density, 

/ = Vf p R n F K O + E U ^ E l + E ^ T > n *> Z ^ + n 'v) + ^K)}> ( 19 ) 
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and the index j runs over all nuclear species and nucleons. Taking j = p, n, we find as usual 
that the Lagrange multipliers, a and /3, are related to the chemical potentials of proton 
(fi p = df /dn p ) and neutron (p, n = df/dn n ) as 

a = /i n , (20) 
= fip-fx n . (21) 

The differentiation with respect to the number densities of nuclei, on the other hand, gives 
the usual relations between the chemical potentials of nuclei \x% and those of proton and 
neutron: 

m = Zi/jp + (Ai - Zi)p, n . (22) 

Differentiating the free energy for nuclei with respective to the number densities of 
nuclei, n,, and employing Eq. fl22|) . we can express as follows: 



/ Zifr + (Ai - Zi)fin -MA 

rii = gi n Qi exp — — — . 23 

V k B T(l - n B /n s ) J 

Note that the only difference from the counter part in the ordinary NSE calculations, 

(ZiUp + [Ai - Zi)p, n - MA 
rii = gi n Qi exp I ^ 1 , (24) 

seems to be the factor (1 — nB/n s ) in the denominator in the exponential. 

This is not true, however, and the calculations for our model are much more involved 
than for the ordinary NSE. In the latter case, all we have to do is to solve the two conservation 
equations, Eq. (TT7|) . for \i p and /x n , since the number densities of all nuclei as well as nucleons 
are expressed by them alone (e.g. Eq. ( T24")) ). In our case, on the other hand, the masses of 
nuclei, Mj, depend on the number densities of nucleons in the vapor, n' p and n' n , which are 
in turn related with \x v and \i n as follows: 

a , - 9f - u' RMF (n> nM + Vn dM ^ n >'n) (25) 

where y^^^ are the chemical potentials of nucleons in the vapor of volume V, which are 
obtained from RMF (see section I2.ip . and the second term that is originated from the 
dependence of the mass of nuclei on the nucleon densities in the vapor and hence is summed 
over all nuclear species is the source of the complication in the solution. As a result, we have 
to solve not only the two conservation equations, Eq. ( ITT)) , but also the relations between 
fip/n and n/ p , n , Eq. ( 125)) . for the four variables, /i p , fi n , n' p and n' n , which then give the number 
density of nuclei, n i? by Eq. (1231) . 
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2.4. Thermodynamical Quantities 



After minimization, we obtain the free energy density together with the abundance of 
various nuclei and free nucleons as a function of pb, T and Y v . As mentioned repeatedly, 
the Coulomb-, surface- and translational-energy contributions to the free energy of nuclei 
are going to zero as the density increases and only the nuclear bulk-energy and free-nucleon 
contributions remain just prior to the transition to uniform nuclear matter. Since both of 
them are entirely based on RMF at this point, the free energy density obtained by the sum of 
them coincides with the one for uniform nuclear matter obtained by RMF. At low densities, 
on the other hand, the excluded-volume effect is negligible and the masses of nuclei agree 
with those of the experimental data. In addition, the RMF free energy for the vapor nucleons 
reduces to the one for the ideal Boltzmann gas of nucleons, since the interactions between 
nucleons become negligible. As a result, the free energy density obtained with the present 
formulation reproduces the well known free energy of the ideal gas composed of nucleons and 
nuclei, the masses of which are given by the experimental data. 

Once the free energy is obtained, other physical quantities are derived by its partial 
differentiations. The baryonic pressure, for example, is obtained by the differentiation with 
respect to the baryonic density as follows: 



where p p ^ F is the contribution from the nucleons in the vapor; p*l ulk , p s ^ r ^ and pf' oul origi- 
nate from the bulk, surface and Coulomb energies of nuclei in the free energy, respectively; 
both pf 1 = riikBT(l — nB/n s ) and pf x = nikBT(n B /n s ) (\og(rii/nQi) — 1) come from the 
translational energy of nuclei in the free energy, with the former being the ordinary pressure 
of nuclei as a Boltzmann gas and the latter arising from the excluded-volume effect approxi- 
mated by the factor (1 — ns/n s ). They are minor at high densities. Note also that p\ ulk and 
p* ur f are normally negligible. The Coulomb-energy contribution, pf' oul , is negative owing to 
the attractive electrostatic interactions among uniformly-distributed electrons and protons 
inside nuclei and is particularly important for ps > 10 n g/cm 3 and Y p > 0.3. If nucleus i is 
in the ordinary droplet phase, it is given as 



Pb = n B [df/dn B ] Ti Ye- f 



(26) 
(27) 




Pi 



Coul 




xi(-«; /3 +ui) 



{■qn e - n' p )(l - Uj) 
(n e - n' p )r]n B 
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whereas in the bubble phase it is expressed as 



Pi 



Coul 





(ryn e - n' p )(l - u<) 
(n e - n' p )r]n B 



(29) 




(30) 



= ^ + E^[{^- g (^)}(l-n B /n s ) 



+{l-n B /n s )Tg' i /g i -dM i /dT , 



(31) 



where g-(T) = [9gj(T)/9T] ye . The first term in Eq. (I3~T|) is the entropy of uniform nuclear 
matter and the second one is the entropy of nuclei as a Boltzmann gas with the excluded- 
volume effect, which is negligible at low densities. The third term accounts for the thermal 
excitations of nuclei with (1 — ns/n^Tg^/ Qi being non-negligible at high temperatures as 
shown later. The last term is originated from the fact that the bulk energy of nuclei is 
temperature-dependent in our formulation and given as follows: 



The contribution of this term is normally negligible except near the nuclear saturation den- 
sity. 

The free energy thus obtained is continuous at any density but not smooth at the 
saturation density. In fact, the pressure, or the density-derivative of the free energy, is 
not continuous at the saturation density. This is most easily understood from Eq. (j27|) . 
The contribution from the excluded volume effect in the translational energy, pf x , in the 
equation does not vanish at the saturation density. The continuity of pressure would have 
been obtained if we had employed a suppression factor (the last term on the right hand 
side of Eq. ( II 5p ) with a higher power. We do not think this is necessary, however, since the 
discontinuity is practically negligible compared with the dominant electron pressure in core- 
collapse simulations. The derivatives of free energy with respect to T and Y p are continuous 
at any density as observed in Eq. (|3T|) for entropy. 



dMi _ dE\ 



(T, n s i, Zi 



(32) 
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Results 



We show the main features of our EOS for selected combinations of densities (ps > 
10 6 g/cm 3 ), temperatures (T = 1, 5, lOMeV) and proton fractions (Y p = 0.1,0 .3, 0.5), empha- 
sizing in particular the differences from the H. Shen's (IH. Shen et al.lll998l ) and Hempel's 
( IHempel et al.ll20101 ) EOS's. The matter compositions are given first. It is also demonstrated 
that the transition to uniform nuclear matter in our EOS occurs in the conventional way, 
that is, either by coalescence of nuclei or by the increase in the density of nucleon vapors. 
Then we turn to the thermodynamical quantities such as free energy, pressure and entropy. 



3.1. Compositions 

The mass fractions of nuclei are shown in the (N, Z) plane with N being the neutron 
number for ps = 10 11 g/cm 3 , T = lMeV and Y p = 0.3 in Fig. [TJ It is clear that nuclei 
are abundant in the vicinities of the neutron magic numbers (N = 28, 50 and 80). In 
fact, the abundance peak arises around (N, Z) ~ (50, 30) in our EOS in this case. For 
the same density, temperature and Y p , on the other hand, the H. Shen's EOS predicts that 
the representative nucleus has (N,Z) = (57.2,33.2). The difference originates mainly from 
the fact that the H. Shen's EOS does not include the shell effect in their Thomas-Fermi 
approximation. At higher densities, heavier nuclei become abundant, which is evident in 
Fig. |2]that shows the mass fractions of nuclei for pB = 10 13 5 g/cm 3 , T = lMeV and Y p = 0.3. 
It should be noted that nuclei with Z > 100 are q uite abundant and, i n fact, are dominant 



in this case. Hence the neglect of these nuclei in IHempel et all (120101 ) cannot be justified 
at these high densities. It should be also mentioned that the nuclear distribution is much 
smoother in this case than that for the lower density. This is because the shell effect is 
reduced in our EOS at high densities, where the nuclear bulk energy is interpolated from the 
one derived experimentally and the one obtained theoretically for uniform nuclear matter, 
with the latter being dominant at this high density. The H. Shen's EOS predicts that the 
most abundant nucleus is the one with (N,Z) = (203.8,89.4), which is still different from 
the prediction of our EOS, (N,Z) = (243,110). These diff erences may ha ve an important 



ramification for the electron capture in the collapsing core ( IHix et all 120031 ). 



In order to expedite comparison with the existing EOS's with SNA, we show in Fig. |3]the 
mass fractions of light {Z{ < 5) and heavy (Zi > 6) nuclei together with free nucleons in the 
vapor for T = lMeV. Irrespective of the values of Y p selected here (Y p = 0.5, 0.3, 0.1), the free 
nucleons dominate other constituents up to pb ~ 10 7 g/cm 3 at this temperature. At larger 
densities (10 7 g/cm 3 < pb < 10 10 g/cm 3 ), some of free nucleons are taken into symmetric light 
nuclei such as deuterons and alpha particles. For densities larger than pb ~ 10 10 g/cm 3 , light 
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nuclei disappear and heavy nuclei are populated dominantly. Near the nuclear saturation 
density, the matter is mostly made of nuclei in the bubble or uniform states. In the case 
of Y p = 0.5, the free protons coexist with symmetric light nuclei and neutron- rich heavy 
nuclei at p# ~ 10 10 g/cm 3 and only symmetric heavy nuclei survive at larger densities. For 
Y p = 0.3, 0.1, on the other hand, most of free protons are contained in nuclei and free neutrons 
coexist with heavy nuclei. In the case of Y p = 0.1, in particular, there is a density regime 
near the nuclear saturation density, where very neutron-rich light nuclei such as the one with 
(Z, A) = (1, 10) exist in the pasta or uniform states as we mention later again. Fig. H] shows 
the matter compositions at a bit higher temperature, T = 5MeV. It is evident for all the 
three values of Y p that free nucleons are dominant up to higher densities, ps ~ 10 12 g/cm 3 . 
It is mentioned incidentally that the mass fractions of light nuclei in our EOS are larger in 
general than those of alpha particles in the H. Shen's EOS. This is partly because the H. 
Shen's EOS has only alpha particles as light nuclei whereas ours takes into account other 
nuclei such as deuterons. The difference somewhat affects the thermodynamical properties. 
For example, the entropy per baryon in our EOS is larger than that in the H. Shen's EOS 
at pb ~ 10 12 g/cm 3 and T > 5MeV owing to larger degrees of freedom in the light nuclei in 
the former as shown in 13.21 

The average mass number, A, of heavy (Z, > 6) nuclei is displayed in Fig. [5j For 
comparison, not only our own results but those of Hempel's and the mass number of the 
representative nuclei, A rep , for the H. Shen's EOS are also shown. It is found that for 
T = lMeV (the left column of the figure), the former two give the average mass numbers 
that are systematically smaller than the H. Shen's EOS at ps < 10 12 g/cm 3 , around which the 
supernova core becomes optically thick and neutrinos are trapped and equilibrated thermally 
and chemically with the matter. It is also clear that the average mass number grows step- 
wise for our and Hempel's EOS's whereas it is monotonic for the H. Shen's EOS. This is 
simply due to the shell effect: nuclei are abundant in the vicinity of the neutron magic 
numbers (N = 28,50,80). At higher densities, ps > 10 12 g/cm 3 , the effect is suppressed in 
our EOS, as mentioned already, because the nuclear bulk energy is interpolated from the 
experimentally derived value and the one obtained theoretically by RMF, which does not 
include the shell effect. Note also that at these densities there appear very heavy and/or 
neutron-rich nuclei, for which no experimental mass data are available and we employ RMF 
for the estimation of their bulk energies (see section l2.2.3p . Heavy nuclei continue to grow up 
until the nuclear saturation density in our EOS, whereas the growth stops at some density 
and nuclei disappear before the saturation density in the Hempel's EOS (see the insets of the 
figure), which is at odds with the conventional idea. This point will be mentioned shortly 
again. 

The discernible difference between our results and the Hempel's at ps ~ 10 12 g/cm 3 for 
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Y p = 0.1,0.3 is attributed to the difference in the masses of neutron-rich nucle i, which we 



calcul ate by RMF whereas iHempel et al.l (120101 ) obtained from the mass table by lGeng et al. 



(120051 ). The dip in the mass number at ps ~ 10 g/cm for Y p = 0.1 in our EOS is artificial 



being produced by the neglect of the shell effect in our mass formula. As nuclei become 
very neutron-rich, there happens a rather rapid shift of abundant nuclei from those with 
experimental data to those without, resulting in the observed change of the mass number. 
The mass number does not increase until higher densities are reached. As shown in the 
right column of Fig. for T = 5MeV, for example, the mass number starts to increase at 
~ 10 13 g/cm 3 . This is because the increase of entropy by the translational degrees of freedom 
of free nucleons is effective to minimize the free energy at this high temperature, whereas 
the decrease of internal energy by the binding energy of nuclei is more advantageous at lower 
temperatures. Another interesting feature at this temperature is the decrease of the mass 
number near the nuclear saturation density for Y p = 0.3, 0.5 even in our EOS. There are two 
reasons for this: one is the formation of bubble nuclei and the reduction of surface energy 
as a result; the other reason is that the less heavy nuclei become, the greater entropy they 
have thanks to their larger translational degrees of freedom and they lower the free energy. 
It should be noted that unlike in the Hempel's EOS, nuclei are not dissociated completely 
in our EOS. 

The rate of coherent scatterings of neutrinos on nuclei is proportional not to the mass 
number of nuclei but to its square as long a s the wave length of neutrino is much larger 



than the nuclear size ([Thompson et al.l 120031 ). As mentioned already, in SNA the rate is 



evaluated by the mass number squared for the representative nucleus, the validity of which 
needs confirmation. In the top panels of Fig. |6]we compare the average mass number squared, 
A 2 , with the square of the average mass number, A 2 , and the mass number squared of the 
representative nuclei in the H. Shen's EOS. In the middle and bottom panels, the standard 
deviations, a a = V A 2 — A 2 , and the dispersion normalized by the average mass number 
squared, a A /A 2 , are also shown. It is intriguing to find that the average mass number 
squared is almost indistinguishable from the square of the average mass number. The reason 
for this small dispersion is that heavy nuclei are rather concentrated in the vicinities of the 
neutron magic numbers. In this sense, SNA is not so bad an approximation in evaluating 
the rate of coherent scatterings of neutrinos on nuclei. It should be noted, however, that the 
mass numbers of the representative nuclei in the H. Shen's EOS are systematically larger 
than the average mass numbers obtained in our EOS. Incidentally, at very high densities 
(Pb > 10 14 g/cm 3 ), most of nuclei are in the pasta phase. 

We demonstrate here that the transition to uniform nuclear matter occurs in one of the 
following two ways depending on the proton fraction of the entire system, Y p . The first case 
corresponds to relatively large proton fractions, Y p > 0.3. At p# > 10 13 g/cm 3 , most of the 
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nuclei in abundance have a similar proton fraction Zi/Ai, which is close to Y p as shown in 
Fig. [2j As the density rises up to the nuclear saturation density, their W-S cells get smaller 
and these abundant nuclei enter the pasta phase almost simultaneously In the bubble phase, 
the volumes occupied by the nuclei in their W-S cells become larger and finally the nuclei 
coalesce with each other. The sequence is schematically illustrated for the nucleus with 
(Z, A) = (6, 20) in Fig. [7J Note that the density in the vapor, or the bubble, remains lower 
than that in the nuclei up to the saturation density in this case. 

The second case occurs for more neutron-rich systems, Y p < 0.3. The proton fractions, 
Zi/Ai, of nuclei having large mass fractions are not necessarily close to Y p of the entire 
system and take various values at ps > 10 13 g/cm 3 . This is due to the fact that more 
symmetric nuclei have greater binding energies. Hence the nuclei are less neutron-rich than 
the vapor in general. As a result of this diversity of the proton fractions of the nuclei in 
abundance, these nuclei reach the pasta phase at different densities. Put another way, nuclei 
of different shapes coexist with substantial mass fractions at a particular density above 
Pb ~ 10 13 g/cm 3 . For example, at Y p = 0.1, ps = 10 13 ' 9 g/cm 3 and T = lMeV, a very heavy 
nucleus with (Z, A) = (48, 228) still remains in the normal phase whereas a very neutron- 
rich light nucleus with (Z, A) = (1, 10) as well as a neutron-rich, slightly heavy nucleus with 
(Z,A) = (6,49) are already in the bubble phase. This light nucleus becomes uniform in 
the W-S cell already at p# = 10 14 ' 1 g/cm 3 whereas the others are still non- uniform. This 
actually corresponds to the formation of a big nucleus with no bubbles inside, which coexists 
with other bubble nuclei. As the density increases further beyond the nuclear saturation 
density up to the point, where the entire system becomes uniform, the density of this big 
nucleus also increases, which we referred to as overcompression in section 12.21 The reason 
for the occurrence of such overcompression is that the binding energies of nuclei that are 
less neutron-rich than the whole system are greater in general and the non-uniform matter is 
energetically favored still at the nuclear saturation density for the whole system. As shown in 
Fig. [HI the transition to uniform nuclear matter for the very neutron-rich light nuclei occurs 
by coalescence with each other just as in the first case mentioned above. On the other hand, 
the transition for the nuclei with (Z, A) = (6, 49) proceeds not by coalescence but by the 
increase in the density of the vapor. Although it is not shown in the figure, it is interesting 
that the very heavy nucleus with (Z, A) = (48, 228) does not reach the bubble phase and is 
in the intermediate pasta phase just prior to the transition to uniform nuclear matter. 



3.2. Thermodynamical Quantities 



Once the minimization is completed, the baryonic contribution to the free energy, to 
which our discussions are limited in this paper, is obtained as a function of ps, T, Y p . All the 
thermodynamical quantities are then derived from it through thermodynamical relations. 
It turns out that most of them are not very different from the ones in the Shen's and 
Hempel's EOS's: for the temperatures, T = 1,5, lOMeV considered in this paper, the three 
EOS's are essentially identical at low densities, ps < 10 8 g/cm 3 , since the matter is simply 
composed of free nucleons as an ideal Boltzmann gas in all the EOS's. At higher densities, 
10 8 ~ Pb ^ 10 13 g/cm 3 , a small difference appears, since some of the free nucleons are taken 
into nuclei, the treatment of which is different among the EOS's. At still greater densities, 
Pb ^ 10 13 g/cm 3 , the difference is more apparent. In the following we will look at the behavior 
more in detail. 

Fig. shows the free energies per baryon as a function of density for the three combi- 
nations of temperature and proton fraction: (T = lMeV, Y p = 0.5), (T = 5MeV, Y p = 0.3), 
(T = lOMeV, Y p = 0.1). In all the cases, the free energies are close to one another among the 
three EOS's. Note that our EOS coincides with the H. Shen's EOS at the highest density as 
it should do. This is also true of the Hempel's EOS although how uniform nuclear matter 
is reached is different between the two EOS's as described in the previous section. Looking 
more closely, we find for T = lMeV and Y p = 0.5 that the Hempel's free energy is larger 
than those of our and H. Shen's EOS's at ps ~ 10 14 g/cm 3 and the hump is followed by a 
steep decline. The difference originates from the fact that the Hempel's EOS includes only 
nuclei with Z < 100 and that unlike our treatment only the Coulomb energy is manipulated 
by using the W-S cell to account for the existence of other nuclei. Their prescription lead to 
the total dissociation of heavy nuclei to free nucleons before reaching the nuclear saturation 
density. The difference between our EOS and the H. Shen's EOS, on the other hand, comes 
from the fact that the latter includes only single representative heavy nuclei. In the case of 
T = 5MeV and Y p = 0.3, our and Hempel's free energies are smaller than that of H. Shen's 
EOS at 10 12 < pb < 10 13 g/cm 3 . Under this condition light nuclei (deuterons in particular) 
are abundant as shown in Fig. HJ The H. Shen's EOS takes into account alpha particles 
alone as light nuclei. The greater degrees of freedom in the light nuclei in our EOS, espe- 
cially the substantial population of deuteron, make the difference observed. The deviation of 
Hempel's free energy from ours at pb > 10 13 g/cm 3 occurs for the same reason as in the case 
of T = lMeV and Y p = 0.5. For T = lOMeV and Y p = 0.1 our free energy is smaller than 
those of the other two EOS's at ps > 10 13 g/cm 3 . The reason is that the existence of pasta 
phase is taken into account in our EOS. In fact, under this condition our EOS predicts the 
co-existence of various neutron-rich nuclei in the pasta phase or in the uniform phase that 
corresponds to a big nucleus with no bubbles inside, and free neutrons (see sectior J3TTl and 
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Fig. [8]). This is in sharp contrast to the matter composed of normal nuclei and free neutrons 
considered in the other two EOS's. 

The pressure is shown as a function of density for three combinations of temperature 
and proton fraction in Fig. [TUJ The three EOS's agree with one another at low densities, 
Pb ^ 10 12 g/cm 3 . The difference appears, however, in the density regime where the pressure 
becomes negative for the H. Shen's and Hempel's EOS's. The baryonic pressure can be 
negative when the Coulomb-energy contribution, which is negative owing to the attractive 
nature of the Coulomb interactions, dominates over the other positive contributions. Note 
that the positive leptonic pressure is much larger than the baryonic one normally and the 
total pressure never becomes negative. For T = 5MeV and Y p = 0.3, our EOS gives positive 
pressures at p B ~ 10 14 g/cm 3 , where the other EOS's predict negative values. This is due 
to the existence of the pasta phase, which our EOS takes into account but the others do 
not. In fact, the absolute values of the Coulomb-energy contribution are reduced when the 
change of the charge distribution in the W-S cell is considered in the pasta phase. As a result 
the pressure in our EOS remains positive. It is also noted that the contribution from the 
excluded- volume effect, "YliiPT ■> gives a non-negligible positive contribution to the pressure. 
At the lower temperature, T = lMeV, even our EOS gives negative baryonic pressures and 
the main difference from the other EOS's lies in the lowest value, p B ~ 10 12 g/cm 3 , of the 
densities that give negative pressures. It is found that our EOS predicts a slightly higher 
density than the others. This is because our EOS gives larger mass fractions of free neutrons, 
which give a positive contribution to the pressure. The difference from the H. Shen's EOS is a 
consequence of the shell effect, which is naturally included in the experimental mass data but 
ignored in the Thomas-Fermi approximation adopted in the H. Shen's EOS. The difference 
from the Hempel's result, on the other hand, comes from the difference in the theoretical 
estimate of the masses of neutron-rich nuclei with no experimental data. In the case of 
T = lOMeV and Y p = 0.5, our and Hempel's pressures are larger than that of H. Shen's EOS 
at p B > 10 13,5 g/cm 3 . This is because the H. Shen's EOS neglects the contribution from the 
translational motions of nuclei, which is not negligible at this high temperature. It is also 
noted that at high densities, ps ~ 10 13 g/cm 3 , the contribution from the excluded-volume 
effect is not negligible, either, the treatment of which is different between our and Hempel's 
EOS's (see section 

The entropy per baryon is displayed as a function of density for three combinations 
of temperature and proton fraction in Fig [Til For high temperatures, T = 5, lOMeV, our 
results are larger than those of H. Shen's EOS at ps ~ 10 12 g/cm 3 owing to larger population 
of light nuclei such as deuteron as well as the inclusion of the entropy from excited states of 
nuclei, gi{T) in Eq. ffl6|) . The difference at p# ~ 10 13 g/cm 3 is originated from the existence 
of heavier nuclei in our EOS than in the others. Thanks to the larger degrees of freedom, 
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the entropy is actually larger in our EOS. At lower temperatures, T = lMeV, the entropy 
per baryon in our EOS is almost equal to those in the other EOS's. 



4. Summary and Discussions 

We have calculated the baryonic equation of state at sub-nuclear densities, which is 
meant for the use in core-collapse simulations and provides the abundance of various nuclei 
up to the atomic number of ~ 1000. The formulation is based on the mass formula and 
the model free energy is constructed so that it could reproduce the ordinary NSE results 
at low densities and make a continuous transition to the EOS for supra-nuclear densities, 
for which we have adopted the H. Shen's EOS in this paper. The experimental mass data 
have been used, whenever available, to obtain the bulk energy of nuclei. For very heavy 
and/or very neutron-rich nuclei with no experimental mass data available, we have utilized 
for the evaluation of the bulk energy the theory used in the EOS for supra-nuclear densities, 
that is, RMF in this paper. At high densities, where the nuclear structure is affected by the 
presence of other nuclei, nucleons and electrons, we have approximated the bulk energy of 
nuclei by interpolation of the value obtained experimentally and the one derived theoretically. 
Assuming the charge neutrality in the W-S cell, we have calculated the Coulomb energy 
of nuclei. Close to the nuclear saturation density, the existence of the pasta phase has 
been taken into account in calculating the surface and Coulomb energies. The free energy 
of the nucleon vapor outside nuclei is calculated by RMF. These prescriptions not only 
ensure the continuous transition to the EOS for supra-nuclear densities wh en the matter 



becom es uniform but also has remedied the peculiar behavior of the EOS by iHempel et al. 



( 120101 ) . in which all nuclei are dissociated completely to nucleons before reaching the nuclear 
saturation density. In our EOS, the transition to uniform nuclear matter proceeds either via 
the coalescence of nuclei or the rise of the density in the nucleon vapor, depending on the 
proton fraction. 

For some representative combinations of density, temperature and proton fraction, we 
have made a comparison of the abundance of nuclei as well as thermodynamical quantities 
obtained in this paper with those provided by the H. Shen's and Hempel's EOS's. It is 
found that the thermodynamical quantities such as the free energy, pressure and entropy 
per baryon, agree well among these EOS's. The matter compositions show noticeable dif- 
ferences, however. The average mass numbers, for example, are systematically smaller at 
Pb ^ 10 12 g/cm 3 in our EOS than the mass numbers of the representative nuclei in the H. 
Shen's EOS. The difference is ascribed to the neglect of the shell effect in the latter EOS. We 
have also demonstrated that very heavy nuclei with Z > 100, which were completely ignored 
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in the Hempel's EOS, are indeed abundant at ps > 10 13 g/cm 3 and lead to the difference in 
the entropy per baryon between the two EOS's. It is interesting that the mass fractions of 
light nuclei such as deuteron cannot be ignored at high temperatures and, in fact, are larger 
than that of alpha particles in the H. Shen's EOS. The average of mass number squared is 
important for the evaluation of the coherent scattering rate and has been frequently approx- 
imated by the mass number squared of the representative nuclei in supernova simulations 
that employ an EOS with SNA. We have shown that the root mean square of mass number 
is very close to the average mass number owing to the shell effect and SNA is justified in 
this sense. It should be noted again, however, that the average masses in our EOS are sys- 
tematically different from the masses of the representative nuclei in the H. Shen's EOS and 
the coherent scattering rate will be reduced in core collapse simulations with our EOS. 

Although we have not made a direct comparison, it is interesting to make some comments 



on ot her EOS's. As mentioned in Introduction, the Lattimer & Swesty's EOS (ILattimer et al. 



199ll ). which is based on the compressible liquid drop model with SNA and a free energy per 



baryon with the Skyrme-type parametrizati on, is another EOS t hat is commonly employed 



in core collapse simulations. Very recently, |G. Shen et al.l (120111 ) published a new EOS for 



use in astrophysical simulations, which employs different theories, the Virial expansion at 
low densities and SNA with the Hartree approximation at high densities. It is noted first 
that since these two EOS's behave quite similarly to the H. Shen's EOS, the results obtained 
for the latter in this paper will be mostly applicable to the former two. For example, 
thermodynamic quantities are not very different although the entropy in our EOS tends to 
be a bit higher than those in others because of the inclusion of various light nuclei in our 
EOS. On the other hand, the step-wise growth of average mass number at low densities is 
common only to our and G. Shen's EOS, since the latter also employs experimental mass 
data. At high densities, the Lattimer & Swesty's and G. Shen's EOS's employ SNA just as 
the H. Shen's EOS. Not including bubble nuclei, the H. Shen's and G. Shen's EOS's predict 
very big droplet nuclei where the Lattimer & Sw esty's and our EOS g ive bubble nuclei. It is 



interesting, however, to mention incidentally that lG. Shen et al.l ( 120 111 ) also found a shell-like 
hollow nuclei, which others did not. 

There is an ample room for improvement in our EOS. The interpolation of the bulk 
energy and the treatment of the pasta phase are entirely phenomenological and need justifi- 
cation or sophistication somehow. Our results show that the shell effect is important for the 
abundance of nuclei. In the present prescription, however, the effect is suppressed maybe 
too strongly at high densities. We are planning to adopt an empirical formula for this effect 
from the literature. With these caveats in mind we are currently constructing an EOS table, 
which will be available for supernova simulations and put in public domain. In so doing, 
electron capture rates, which are consistent with the abundance given by the EOS, should 
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be also included. In our formulation, it is quite simple to change the EOS for supra-nuclear 
densities. In fact, all we need is the free energy density of uniform nuclear matter predicted 
by the particular EOS. The combination with another EOS is under progress at present. 
Different RMF parametrizations will be also discussed in our forthcoming paper. 

The authors gratefully acknowledge the contribution of A. Ohnishi and C. Ishizuka 
to the core idea. We would like to thank M. Takano for useful information on the mass 
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Education, Culture, Sports, Science and Technology (19104006, 21540281, 22540296) and 
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neutron number N 

Fig. 1. — The mass fractions of nuclei in the (N, Z) plane for Pb = 10 11 g/cm 3 , T = lMeV 
and Y p = 0.3. The cross indicates the representative nucleus for the H. Shen's EOS under 
the same condition. 
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Fig. 2. — The mass fractions of nuclei in the (N, Z) plane for p# = 10 13,5 g/cm 3 , T = lMeV 
and Y p = 0.3. The cross indicates the representative nucleus for the H. Shen's EOS under 
the same condition. 
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Fig. 3. — The mass fractions of free protons (dotted lines), free neutrons (dashed lines), 
light nuclei (Z < 5, long dashed lines) and heavy nuclei [Z > 6, solid lines) as a function of 
density for T = lMeV and Y p = 0.5, 0.3, 0.1 from top to bottom. 
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Fig. 4. — The mass fractions of free protons (dotted lines), free neutrons (dashed lines), 
light nuclei (Z < 5, long dashed lines) and heavy nuclei [Z > 6, solid lines) as a function of 
density for T = 5MeV and Y p = 0.5, 0.3, 0.1 from top to bottom. 




Fig. 5. — The average mass number, A, of heavy nuclei with Z > 6 for our EOS (solid 
red lines) and the Hempel's EOS (dotted green lines) together with the mass number of 
representative nucleus for the H. Shen's EOS (dashed blue liens) as a function of density for 
T = lMeV and Y p = 0.1, 0.3, 0.5. The insets are the close-ups of the high density regimes. 



- 31 - 




lo Kio P,le/™ 5 ] log, pjg/cm 3 ] log,,, pjg/cm 3 l 





Fig. 6. — The square of mass numbers (top), the standard deviation of mass number, 
a a = \/ A 2 — A 2 , (middle) and the dispersion normalized by the average mass num- 
ber squared, cr\/A 2 , (bottom) of heavy nuclei with Z > 6 for T = lMeV and Y p = 
0.1 (left), 0.3 (middle) and 0.5 (right). In the top panels, the solid and dotted lines show the 
average mass number squared, A 2 , and the square of average mass number, A 2 , in our EOS, 
respectively, whereas the dashed lines display the mass number squared of the representative 
nucleus for the H. Shen's EOS. 
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Fig. 7. — Schematic pictures of the number densities of proton (dashed red lines) and neutron 
(solid blue lines) in the W-S cells for the nucleus with (Z, A) = (6, 20) at T = lOMeV, 
Y p = 0.3 and ps = 13.5, 14.0, 14.1g/cm 3 from top to bottom. The dotted black lines indicate 
the cell boundaries. 
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Fig. 8. — Schematic pictures of the number densities of proton (dashed red lines) and neutron 
(solid blue lines) in the W-S cells for the nuclei with (Z, A) = (6, 49) (left) and (Z, A) = 
(1, 10) (right) at T = lMeV, Y p = 0.1 and p B = 13.5(left), 13.8(right), 14.0, 14.1g/cm 3 from 
top to bottom. The dotted black lines indicate the cell boundaries. 
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Fig. 9. — The free energy per baryon for our (solid lines), H. Shen's (dashed liens) and 
Hempel's (dotted lines) EOS's as a function of density for (T = lMeV, Y p = 0.5), (T = 
5MeV, Y p = 0.3) and (T = lOMeV, Y p = 0.1) from top to bottom. 
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Fig. 10. — The baryonic pressure for our (solid lines), H. Shen's (dashed liens) and Hempel's 
(dotted lines) EOS's as a function of density for (T = 5MeV, Y p = 0.3), (T = lMeV, 
Y p = 0.3) and (T = lOMeV, Y p = 0.5) from top to bottom. 
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Fig. 11. — The entropy per baryon for our (solid lines), H. Shen's (dashed liens) and Hempel's 
(dotted lines) EOS's as a function of density for (T = lOMeV, Y p = 0.5), (T = 5MeV, 
Y p = 0.3) and (T = lMeV, Y p = 0.3) from top to bottom. 



